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1. Introduction 



There has been considerable interest in the idea that our universe may be a 3-brane 
embedded in a higher dimensional bulk space. This possibility has been motivated 
by ideas in string/M-theory. All the usual standard model fields are confined to the 
brane, while gravity can propagate in the bulk. So far, gravity experiments have not 
detected any higher dimensional effects, so if this 'brane world' scenario is to be a 
model of our universe, they must be suppressed. In a scenario considered by Randall 
and Sundrum a positive tension 3-brane is embedded in a five-dimensional bulk 
space with a negative cosmo logical constant. If the brane tension and bulk energy 
density are chosen appropriately, the standard four-dimensional gravity will arise 
from graviton bound states on the brane. The gravitons which are free to propagate 
in the bulk give corrections to this. 

This idea was extended to an evolving universe by Binetruy, Deffayet, Ellwanger 
and Langlois |2|, [|. As well as the usual brane energy density there is also a brane 
cosmo logical constant. If this is chosen appropriately it will cancel the bulk energy's 
contribution on the brane and give a cosmology whose late time Friedmann equation 
is approximately the same as in the standard Friedmann- Robert son- Walker universe. 
Thus the brane cosmology predictions for the current evolution of the universe are 
as normal. However, the initial evolution of the universe is dominated by the non- 
standard terms arising from brane effects. 

In the Randall-Sundrum brane world scenario the bulk contains only a cosmolog- 
ical constant. String/M-theory suggests that it will also contain scalar fields, which 
are free to propagate through the bulk ||. The evolution of these fields may have 
interesting cosmological effects. For example, if the bulk energy density varies with 
time, the effective brane cosmological constant will also be time-dependent. This 
could give alternative mechanisms for inflation and quintessence. 

Solving the field equations for a brane world cosmology with scalar fields is 
far more complicated than in the constant bulk case. One approach is to express 
the five-dimensional Einstein equations in terms of four-dimensional tensors on the 
brane ||, and then use boundary conditions to determine them ||. Unfortunately 
this approach does not fully determine all the tensors. Progress can be made by 
making assumptions about the bulk solution 0, 0, §, but it is not obvious that 
such assumptions are justified, or self-consistent. To determine which solutions of 
the four- dimensional Einstein equations are allowed, the full five-dimensional theory 
must be solved. 

Full bulk solutions for static brane worlds have been found in some special cases. 
If the scalar potential takes a supergravity-like form, the field equations can be 
reduced to first order equations, which can be solved with relative ease || ID| . In 
this paper I will combine that method with the work of Binetruy et al. to give 
cosmological solutions. 
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Conventions and background are outlined in section |2]. In section I extend 
the results of ref. to a special class of brane world cosmologies with bulk scalar 
fields. In section [| the supergravity inspired static solution generating technique 
is generalised to these brane cosmologies. Some examples of this type of solution 
are considered in section |5], and their bulk behaviour is discussed in section [)} The 
results are summarised in section [7]. 



2. Five-dimensional Brane Worlds 

We will consider a theory with the effective five-dimensional action 

s = I d*x^-g\^- 2 R- l -(y<p) 2 -V {<(>))- I d^V^h{- 2 [K] + C h {4>)\ . 

./bulk I ZK Z 2 J Jbrane I K z J 

(2.1) 

qab is the bulk metric, and Jiab is the induced metric on the brane. They are related 
by Hab = Qab — NaNb, where N A is the unit vector normal to the 3-brane. The 
covariant derivatives corresponding to qab and Kab are Va and Da respectively. The 
Gibbons-Hawking boundary term, [K], is required to give a well-formulated action. 
[K] is the change in K = H ab Kab across the brane, and Kab = h^hgV^c^D) is the 
extrinsic curvature, k 2 = 87r/M|, where M5 is fundamental five-dimensional Planck 
mass. £b{4>) is the Lagrangian density for all the fields on the brane. It includes 
the 'ordinary' standard model fields which are confined to the brane, and the brane 
potential energy for the bulk scalar 0. 

Varying the action gives the bulk equations of motion 



dV 



(2.2) 



G A B = k 2 T a = k 2 (V^0V B - ^[|(V0) 2 + Vm) . (2.3) 
The brane part of the action gives the jump conditions 0] 

NAV M = ' ^ 

L J 0(p 

[K AB - Kh AB ] = -kH AB , (2.5) 

where tAB is the energy momentum tensor corresponding to Cb(<p). 
The Codacci equation relates a projection of Gab to Kab 

N A G AB hE = D A K A - D C K . (2.6) 

Combining this with (|2.3|) and substituting it into the jump conditions (|2.4j) and 



implies 



D a4 = - 5 4t D c<P ■ (2-7) 
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This is the effective energy-momentum conservation equation on the brane. 
Following ref. [§] we consider a five-dimensional metric of the form 



ds 2 = g AB dx A dx B = -n 2 (t, y)dt 2 + a z (t, y)Vt ij dx'' dx J + b 2 (t, y)dy 2 (2.8) 

where is a three-dimensional metric of space with constant curvature k = — 1, 0, 1. 
The coordinates along the brane are denoted x^ and y = x 5 is the fifth bulk coordi- 
nate. The brane is located at y — 0. Differentiation with respect to y and t will be 
denoted by primes and dots respectively. The brane's normal is N A = S^b -1 , and 
the extrinsic curvature is = diag(n'/(no), <5*a'/(a&), 0). 

Assuming that the brane matter is describable by a perfect fluid, 
t% = diag(— ph,Pb,Pb,Pb)- The jump conditions (|2.4j ) and (|E5| ) are then 



(2.9) 



Pb , (2.10) 



a b 3 

(3p b + 2p b ), (2.11) 



\n'\ n 2 



a b 3 

The subscript denotes quantities at y = 0, i.e. on the brane, and [Q] = Q(y = + ) — 
Q(y = 0") is the change in Q across the brane. If we assume that the bulk solution 
is symmetric across the brane, i.e. invariant under y «-> — y, then Q(0~) = — Q(0 + ) 
and ( |2.9H2.11| ) can be used to determine cp', a' and n' at y — 0. 
The energy-momentum conservation equation (|2.7|) becomes 



Pb + 3^(p b +p b ) = ^0 o . (2.12) 

The right-hand side of this equation gives the amount of energy non- conservation of 
the matter fields on the brane. 



3. Brane Cosmology with Scalar Fields 

Using the Gauss-Codacci formalism it is possible to express the five-dimensional 
Einstein equations in terms of the effective four- dimensional Einstein tensor, the 
extrinsic curvature, and the electric part of the Weyl tensor p. Using the jump 
conditions the extrinsic curvature is found in terms of the brane energy-momentum. 
The Weyl tensor cannot be found in this way, so another approach is needed to find 
full bulk solutions. I will consider a special class of solutions for which the five- 
dimensional field equations can be simplified and solved using generalisations of the 



methods of refs. 101 and ITlOl . 
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In order to simplify the Einstein equations, it is useful to introduce the quantity 

Ht , v) = -(±)\(^) . (3.1) 
V an J \ab ) 

In the bulk the components of the Einstein tensor are related to T by 

G°o ~ ~,G\ = 7^d y (a^) - 3 4 (3.2) 
a' 2a A a' y a 2 

G\ - -G\ = -^rd t (a^) - 3 4 ■ (3.3) 
a 2a a a a z 

This extends a simplification used in ref. ||. In that paper the bulk energy was taken 

to be constant. The left hand sides of ( |3.2j ) and (!0) are then equal, and it is simple 

to show that they imply T = J 7 (a). The Einstein equations can then be solved with 

relative ease. 

When scalar fields are introduced it is clear, after substituting in ( |2.3|) , that the 
left hand sides of ( |3.2| ) and (|3.3| ) are not generally equal. However if we consider a 
special class of solutions with = 0(a), the two expressions will be the same. Hence, 
as with the constant bulk case, T' ja! = T '/a, and so T = JF(a) too. Both ( |3.2|) and 
^3f) then reduce to 

dT _ ( dd>\ 2 k 2 , , 2k , 

The other diagonal components of the Einstein equation (Gij = k 2 ^) are automat- 
ically satisfied as a result of the Bianchi identities. 

The standard brane cosmology discussed in ref. is a simple example of this 
type of solution. It has = const., V = const. < and 

6 a 2 a 4 

where C is an arbitrary integration constant. 

On the brane the effective four-dimensional Friedmann equation can be deter- 
mined from T . Assuming the Z 2 symmetry, ( |3.1| ) and ( |2.10|) imply 

The time has been scaled so that n = 1. t is thus the proper cosmo logical time on 
the brane. 

When = 0(a), the field equation (|2.2j ) reduces to 



V = a T -— + - — 

V da z a da , 



3 a 



T = %- < 3 ' 7) 

da do 
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Substituting Gab = k 2 Tab and ( |3.4j ) into 



T a 



da 



+ 



dT_ 

da 



+ AT + — T a 



da 



gives 

2" 

a 



dV 



da 



. 



(3.f 



Thus the original partial differential field equations have been reduced to two ordinary 
differential equations. 



4. Supergravity-style Solutions 

For certain choices of V it is possible to generate solutions to the field equations 
using simpler first order equations. If k = and V is a supergravity-style potential 



v{4>) = o 



1 fdW 



-W 2 
6 



(4.1) 



then the field equations 



and 



are satisfied if 



K 



3 dW 



(4.2) 



(4.3) 



and 10 , where 



da k 2 W d(f) 
These equations are a generalisation of those discussed in refs. 
static (if: independent) brane worlds were considered. The corresponding first order 
equations discussed in those papers can be obtained by setting b — 1, n — a — e A ( y \ 

If we were dealing with a supergravity theory then W ((f)) would be the super- 
potential. However supergravity is not required for solutions of ( |4.1| - |4.3| ) to satisfy 
fl3.4p and ( |3.8| ). All that is needed is the above form of the potential. For conve- 
nience I will refer to W as the superpotential whether the underlying theory includes 
supergravity or not. 

The Friedmann equation ( p.6|) becomes 



al 



K 

36 



[P 



(4.4) 



Using ( f4.3| ) the rate of change of on the brane can also be found. 




(4.5) 



The equation ( [4.3[ ) relates the two jump conditions Q2.9[ ) and ( |2.10| ). If the bulk 
is Z2 symmetric, the jump conditions are only consistent if 



5C h 




(4.6) 
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Thus the energy on the brane is not conserved in this type of brane cosmology In 
order for Ch is to satisfy the above condition, there must be dependent couplings 
in the standard model Lagrangian. Some of the models discussed in ref. || (or 
generalisations of them) may be suitable. 

Combining ( |4.6| ) with ( |2.12| ) gives the new energy conservation equation 



Pb + 3— (p b + Pb) = TTTpb • (4.7) 
a Wo 

Thus if Wq/Wq is negative, energy will leak off the brane. If it is positive, energy 
will be pushed out of the bulk and onto the brane. 

In the usual brane cosmology, the energy density pb consists of a cosmological 
constant part and an ordinary matter part. The simplest generalisation of this is 

p b = \W + W p , (4.8) 

where p is proportional to the energy density of the ordinary matter. It has the 
same evolution as in the standard cosmology, i.e. p oc a~ 3 ( 1+w > if p = wp. More 
complicated solutions of (|4.7f) will also exist, but I will not consider them here. 



The effective Friedmann equation (4.4) then becomes 



Qj(\ AC 



1 = ^ (A 2 - 1) + p 2 + 2Ap , (4.9) 
% 3o L J 

It resembles the usual brane cosmology, but with a time dependent gravitational 
coupling strength. If W tends to a suitable constant at large t then a conventional 
brane world Friedmann equation will be obtained. If A = 1 this will give the same 
late time evolution as the standard cosmology, i.e. cio/ a o ^ P- 

By considering the brane Friedmann equation ( |4.9| ), we can rule out many choices 
of W. If W does not tend to a constant, then the solution's cosmological evolution 
will not generally be compatible with observations. It is still possible that such solu- 
tions will agree with the standard cosmology at some intermediate time. However, we 
expect the typical time-scale for these solutions to be of order Mg -1 , so without signif- 
icant fine-tuning of the solution the period of acceptable evolution will be too short. 
Even if the Friedmann equation does have appropriate late-time behaviour, there is 
no guarantee that the effective gravitational interactions of particles on the brane 
will take the standard form. By considering the effective four-dimensional gravity, 
we will be able to further constrain W. However in this paper I will concentrate on 
the implications of the non-standard Friedmann equation. 

It is worth noting that (|4.1|-fO|) still satisfy equations ( |3.4j ) and fl3.8| ) if we change 



the signs of T and V. In this case (ft.4p becomes 



J = sW + rfl. MO) 
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This can also be rewritten as ( f4.9| ), but with (A 2 — 1) replaced with (A 2 + 1). Thus 
these alternative solutions will not give the correct late time behaviour, because of 
the large effective cosmological constant term. 



5. Solutions for Exponential Superpotentials 



As an illustration I will consider superpotentials which are the sum of two expo- 
nentials. These are well motivated by string/M-theory. A time-independent ex- 
ample of one of them is discussed in ref. ||. With the aid of the transformation 
4> — > ±0 + const, they can all be reduced to 



W = c 

with s = ±1 and ai > \(X2\, or if a 2 

W = 



-ai<i 



3 O2 



a i 



+ 8- 



Q-2 



a i 



s 

+ - 

K 



(5.1) 



(5.2) 



Since the above superpotentials are loosely inspired by string theory, it is natural for 
the coefficients «j to be of order k, as in ref. ||. 

If s = +1 and «2 ^ then W has a minimum at <fi = 0. If sa 2 < then W has a 
zero at <p = cp* = (cti + a^) -1 ln(— sa 2 /«i), or cp* = —a^ 1 lncti if s — — 1 and ct2 = 0. 

The corresponding potentials obtained from fl4.1|) are 



v = c - 



4k^ 
3a? 



-2ai0 



1 



4k 2 
3a 9 2 



e 2a 2 4> _ 2s 1 + 



Ak 2 



3ai«2 



3 (a 2 -ai)</> 



and (for a 2 = 0) 



1 - 



4k 2 



-2ai<; 



-2s- 



-ai<, 



, (5.3) 



(5.4) 



3af J 3ai 

If V is to be bounded from below then only some values of the parameters are allowed. 
The solutions of ( ^4.3| ) when W is given by ( |5~T|) are 

,..2 



ln(a/a* 



In 



3aia; 2 

where a* is an arbitrary constant. 

When a 2 = 0, we will use the superpotential 
solved by 



e ai * - se~ a2q 



(5.5) 



!2l) instead. In this case 



is 



ln(a/a* 



/v 



3«i 



[K(t> + se 



(5.6) 



The evolution of each of these solutions falls in to one of six cases, one for each choice 
of s and the sign of a 2 - Since we are trying to produce a big bang style cosmology, 
we will be interested in solutions which start at a = 0. There are also some solutions 
with = const, and dW/dip = 0, but these are effectively the same as the constant 
bulk energy case ( |3.5| ), so I will not consider them further. 
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5.1 Superpotential with a minimum, ct 2 > 0, s = +1 

The expression (|5.5|) gives two solutions in this case. As a goes from to oo, rolls 
down from either +oo or — oo to the minimum of the superpotential at = 0. While 
it is rolling down, W/W < 0, and so (f4.7|) implies that energy leaks off the brane. 

The approximate behaviour of and T in terms of a can be found when a is 
small or large by using series and asymptotic expansions. Initially, when a « 0, 

0^ sgn(0)^ln(a*/a) , (5.7) 
36af \a J 

where the index % is equal to 1 if < and 2 if > 0. Substituting these expressions 
into ( [4.9| ), we find that a oc 

t 1 /(4+3a?/K 2 )_ initia j 

expansion is thus faster than in 

the usual brane cosmology. 
Later, when a is large, 

1 /a*\ 3aia2 / K2 , , 

(5.9) 



Oil + Q!2 V a 



^c 2 



(ai + a 2 ) + aia 2 — 
a 



Thus when t (and hence a) is large, VT 2 — > const, and so the standard cosmology will 
be obtained from (|4.9| ), with order paQ 01 " 2 ^ corrections to the Friedmann equation. 
Since W — > at late time, the amount of energy leakage from the brane also tends 
to zero, as is required to agree with the standard cosmology. 

5.2 Superpotential with a zero, oli > 0, s = — 1 

In this case rolls from either +oo to — oo, or vice versa. At both ends of the solution 
a — » 0. We see from ( |5.5| ) that a is bounded. 

Initially W/W is negative, so energy drains off the brane [see ( [4.71) 1. This contin- 
ues until = 0*, where W = 0. At this point there is nothing left on the brane. This 
is also the point at which the scale factor a reaches its maximum value (a max ). After 
this the universe re-collapses. W/W is now positive, so energy is sucked onto the 
brane. This continues until |0| and \ W\ reach infinity, at which point a has returned 
to 0. 



The approximate behaviours of and J 7 near a = are given by (|5.7| ) and (|5.8| ). 
Near a = a max 

k 2 

ln(a/a max ) w — - (0 - 0*) 2 , (5.11) 
o 

K 2 c 2 e~ 2ai ^* 

T w — -5 (ai + a 2 ) 2 In (a max /a) . (5.12) 



6«i 



S 



By substituting ( |5.8| ) and ( |5.12| ) into (|4.9|) , a finite upper bound can found for the 
time taken to reach a max . While a = there, (|4.5|) shows that ^ and so 
continues to roll down (or up) the superpotential. Thus the universe expands to its 
maximum size and then re-collapses, all in finite time. For «j ~ k, the time taken 
for the universe to re-collapse will be of order M^ 1 . This very short lifetime for the 
universe, and the way in which all energy (and hence matter) flows off the brane, 
mean that this model is completely incompatible with the standard cosmology. 

The re-collapse of the universe is a common feature of solutions which pass 
through zeros of W. If W changes sign and dW/dip ^ 0, then (|4.3|) implies that 
d(p /da changes sign at this point. Equation (|4.5|) ensures that O does not, provided 
p\ > Wq . Thus d must also change sign. This is consistent with the expression for 
do, (fOl). Once do starts to decrease, (|4.4|) implies that it will continue to do so. 



5.3 Positive monotonic superpotential, < 0, s = — 1 

There is just one solution for these parameters. The scalar field starts at — oo when 
a = 0, and then rolls down to +oo as a — > oo. The approximate large and small a 
behaviour of and T are given by ( |5.7|) and (|5.8|) . While the initial behaviour of 
this solution is similar to that of section |5.1] , T — > at late time. Thus the effective 
gravitational coupling tends to zero as t — > oo. 

This type of superpotential has been suggested as a possible candidate for 
quintessence ||. If we set p = and A 2 ^ 1 in ( |4.9| ), we get a solution with a cosmo- 
logical constant which decreases with time. Unfortunately when the matter fields are 
re-introduced, the solution is incompatible with the standard cosmology, due to the 
asymptotically vanishing gravitational coupling. Thus, this type of superpotential is 
not suitable after all, at least not for the cosmological generalisation of supergravity 
solutions discussed in this paper. 



5.4 Superpotential with minimum and zero, «2 < 0, s = +1 

For these parameters a = when = and = — oo. This gives a total of three 
solutions. can roll from — oo to 0, or vice versa. In both cases a —>■ at the 



start and end of the solutions. These are similar to the solutions of section |5^, in 
that the universe expands to a maximum value and then re-collapses. As before the 
maximum value is reached when W = 0. 

The other solution, with > 0, starts at = and a = 0. As a increases to oo, 
tends to +oo. The approximate behaviour of and T there is given by (|5.9|) and 



( [5.1 0| ) near = 0, and is asymptotically given by ( |5.7| ) and ( |5.8| ) when |0| is large 



Thus the > solution has W 2 — > at large t (and a), which is not cosmologically 
acceptable. 
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5.5 Special cases a 2 = 0, s — ±1 



Asymptotic and series expansions of 
large and small a 

3a± 



K 



1|) give the approximate forms of and T for 
In (a/ a*) < 



3«is 



(5.13) 



ln(a/a*) 



> 



K 2 C 2 

36 



3af//c 2 



n 2 



+ S 



+ 1 



< 



> 



(5.14) 



3a? ln(a/a*) 

Both these cases have rolling down the superpotential from — oo to +oo. If s = — 1 
the system will pass through a zero of the superpotential, and as with the solution 



in section 5.2, a reaches its maximum value there, and then starts to decrease. 



The s = +1 case resembles the solution in section 5.1, in that T — > const, at late 



t (large a). Thus this type of solution also provides a potentially viable cosmological 
model. The initial expansion of the universe is a oc 

£l/(4+3a?/ K 2 )_ The late time 

Friedmann equation resembles that of the usual brane cosmology with order pj In a 
corrections. 

5.6 Summary 

The solutions considered in this section can be divided into two distinct categories. 
Firstly there are the solutions which pass through a zero of the superpotential. In this 
case the universe expands until W = 0. Up until this point energy is lost from the 
brane into the bulk. By the time W = there is nothing left on the brane. After this 
the universe contracts, and energy flows back onto the brane. This continues until 
the universe re-collapses to a = 0. This scenario appears to have little resemblance 
to the standard cosmology. The solutions described in sec. [572], sec. |57J] (first two 
cases), and sec. |5T 



(s = —1) fall into this category. 
Of course it is still possible that these solutions agree with the standard cosmol- 
ogy at some intermediate period of their evolution. However, unless the parameters 
are fine-tuned so that cti/ k are extremely small, this period will be far to short to 
agree with observation. 

From a cosmological point of view, the other solutions are more promising. In 
this case <fi rolls down to a minimum of |W|, without passing through W — 0. As 
the minimum is approached, a — ► oo. Energy also leaks off the brane in this case, 
but the rate approaches zero as a —>■ oo. If the final value of W is finite, then the 



standard cosmology can be obtained late in the evolution. The solutions in sec. O 
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and sec. |T5] for s = +1 fall into this category, thus the only viable cosmo logical 
solutions with a superpotential of the form (|5.1|) or (|5.2j ) are those with s = +1 and 
«2 > 0. The solutions of sec. |0| and sec. ^]4| (third case) also have rolling down to 
a minimum, but since it has W = there, the energy remaining on the brane tends 
to zero as a — > oo. 

It should be noted that although the solutions with a minimum of \W\ give the 
correct Friedmann equation (and hence agree with many observations), their other 
properties could still conflict with experimental results, and so they may not be 
acceptable after all. For example, they may give the wrong effective gravitational 
forces for objects on the brane. 

6. Bulk solutions 

So far I have not determined the y dependence of the metric (or <fi). Following 
ref. 0j, this can be found by solving the non-zero off-diagonal component of the 
Einstein equations. When = 0(a) this is 



n'a ba' a' 



,2 



G 05 = 3 ( — + f- - - ] = k 2 T 05 = K 2 a'd ^ . (6.1) 
na ba a \ da 



We can assume that b = 1, in which case Q6.1|) is solved by 



a _ , . h} k 2 



where the function {3 is independent of y. If V has a supergravity-like form (|4.1| ), 
the above expression simplifies to 

- = /9(t)^W(0) • (6.3) 
n 6 



Since we want n = 1, fl4.4Q implies that /3 = aoy pI/Wq — 1. 
Substituting (|6.3|) into the definition of JF (|3.1|) 



(a) 2 = ^(/5 2 + « 2 )^ 2 (0). (6.4) 
36 

a must also satisfy a = a at y = and the jump condition (|2.10|) . Since we are 
interested in a Z 2 symmetric solution, a will be a function of \y\. 

The alternative solutions discussed at the end of section [| satisfy 

(a') 2 = ^(P 2 - a 2 )W\<P) (6.5) 
36 

instead of (|6.4|). Solutions of this equation will have a bounded in the bulk. Since 
they do not have the correct late time cosmological evolution, I will not consider 
them further. 



11 



The corresponding time-independent bulk solutions can be obtained by setting 
Pb = Wo, so that (3 = 0. Finding analytic solutions of ( |6.4|) is not generally possible, 
although solutions for a few special cases can be found. 

If we take a\ = a 2 = k/v3 then the superpotential used in section jO] simplifies 

to 



W 



2c 



a i 



cosh(ai0) 



zy/3 



+ 4 



1/2 



The bulk behaviour of a is 
1 



K2pb sinh(/,M)+^ Pl 



3/i 



2Wi 



(cosh(/iy) 



18/i 2 



(6.6) 



(6.7) 



where p = 2ck/\/3. a 2 has zeros if pb > W^o > 3p/ft 2 , which is true for cosmo logical 
solutions. The position of these singularities varies with time. Unfortunately <ft', <fi 
and W are all singular at these points, and so Tab is too. Thus the singularities 
at a = are naked curvature singularities. This problem could be resolved if we 
have a compact bulk, and one or more other branes at suitable distances from our 
brane TT\. This separation of the branes will then be time dependent. In order for 
this scenario to work, the other brane(s) must have negative energy densities, which 
may lead to instabilities. 

The bulk behaviour for the model discussed in section |5.2| can also be found 
when «! = a 2 = k/\/3. This time the bulk solution of a 2 is given by fl6.7Q , but with 
p replaced by ifi. As before there will be bulk singularities. 

Finally, for ol\ = = «/\/3, the bulk behaviour of the solution in section [5]3 

is 



1 



k 2 k 4 N 

Y P b \y\ + r w 2 y 2 



(6.S 



Again a 2 has zeros in the bulk when pb > Wo, and since Tab is also divergent there, 
they are curvature singularities. In this case W oc 1 /a so the singularities occur when 
\W\ — ► oo. This is also the case with the other two examples above. This suggests 
that solutions which have W bounded may be curvature singularity free. 



7. Conclusions 

By considering solutions of the special form = 0(a) it is possible to significantly 
simplify the field equations for a brane world cosmology with a scalar field. Further- 
more, if the theory has a supergravity style potential, the simplified equations can be 
further reduced to three first order differential equations. This gives a generalisation 
of time-independent supergravity brane worlds to brane cosmologies. For simplicity 
I have considered solutions with Z 2 symmetric bulks. 

Investigation of examples of these solutions suggests that the evolution of these 
brane cosmologies starts and ends at points with either dW/d<p = or W = ±oo (W 
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being the superpotential). If the solution passes through a point with W = 0, the 
universe will stop expanding, and re-collapse. The only solutions which avoid this 
re-collapse are those where <fi rolls down to a minimum of \W\. 

In order for these Z 2 symmetric, supergravity solutions to be self-consistent, the 
matter fields on the brane must couple to the bulk scalar in such a way that energy 
is not conserved on the brane. In the case of solutions which approach a minimum of 
\W\, the amount of leakage tends to zero, and the standard cosmology is obtained at 
late times. With the other solutions, all the energy will have leaked off by the time 
W = is reached. After this energy flows back onto the brane until it re-collapses. 

The bulk behaviour of the metric (and 0) can also be determined. Unfortunately 
for the exponential superpotentials considered in this paper, they all have naked 
curvature singularities in the bulk. The problems occur at points where W = ±oo. 
This suggests that a solution which only passes through finite values of W may 
not have naked bulk singularities. The singularity problem could also be resolved 
by adding additional branes to the bulk. Like the singularities, the position of the 
branes would vary with time. 
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